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Abstract 

The reactions e"'"e~ — >■ e"'"e~AX and e"'"e~ e"'"e~AX are studied using data 
collected at LEP with the L3 detector at centre-of-mass energies between 189 and 
209 GeV. Inclusive differential cross sections are measured as a function of the 
lambda transverse momentum, pt, and pseudo-rapidity, rj, in the ranges 0.4 GeV < 
Pt < 2.5 GeV and \ri\ < 1.2. The data are compared to Monte Carlo predictions. 
The differential cross section as a function of pt is well described by an exponential 
of the form Aexp{—pt/ {pt))- 
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1 Introduction 



Two-photon collisions are the main source of hadron production in high-energy e+e^ inter- 
actions at LEP via the process e+e" —>■ e'^e~j*j* — > e^e~hadrons. The outgoing electron 
and positron carry almost the full beam energy and their transverse momenta are usually so 
small that they escape, undetected, along the beam pipe. In this process, the negative four- 
momentum squared of the photons, Q^, has an average value (Q^) — 0.2 GeV^ and they can 
therefore be considered as quasi-real. In the Vector Dominance Model (VDM), each photon can 
transform into a vector meson with the same quantum numbers, thus initiating a strong interac- 
tion process with characteristics similar to hadron-hadron interactions. This process dominates 
in the "soft" interaction region, where hadrons are produced with a low transverse momentum, 
Pf Hadrons with high pt may also be produced by the QED process 77 qq (direct process) 
or by QCD processes originating from the partonic content of the photon (resolved processes). 

The L3 Collaboration has previously studied inclusive tt". Kg [1] and charged hadron [2] 
production in quasi-real two-photon collisions. In this Letter, the inclusive A and A production^^ 
from quasi-real photons is studied for a centre-of-mass energy of the two interacting photons, 
VF-y-y, greater than 5 GeV. The results are expressed in bins of transverse momentum, p^, and 
pseudo-rapidity, r], in the ranges 0.4 GeV < pt < 2.5 GeV and \ri\ < 1.2. The rj range is further 
divided in two subsamples with 0.4 GeV < < 1 GeV and 1 GeV < pt < 2.5 GeV. The data 
sample corresponds to a total integrated luminosity of 610 pb~^ collected with the L3 detector 
[3], at e"^e~ centre-of-mass energies i/i = 189 — 209 GeV, with a luminosity weighted average 
value of ^/s = 198 GeV . Results on inclusive A production for a smaller data sample at lower a/s 
were previously reported by the TOPAZ Collaboration in the range 0.75 GeV < pt < 2.75 GeV 
[4]. The HI Collaboration investigated the A photoproduction at HERA [5]. 

2 Monte Carlo simulation 

The process e+e" e+e' hadrons is modelled with the PYTHIA [6] and PHOJET [7] event 
generators with two times more luminosity than the data. Both generators include the VDM, 
direct and resolved processes. In PYTHIA, including these processes for each photon leads 
to six classes of interactions. A smooth transition between these classes is then obtained by 
introducing a transverse momentum scale. The two-photon luminosity function is based on the 
EPA approximation [8] with a cut at the mass of the rho meson. The SaS-lD parametrisation 
is used for the photon structure [9]. PHOJET rehes on the Dual Parton Model [10], with a soft 
and a hard component. The two-photon luminosity functions is calculated in the formalism of 
Reference 8. The leading-order GRV parametrisation is used for the photon structure [11]. For 
both programs, matrix elements are calculated at the leading order and higher-order terms are 
simulated by parton shower in the leading-log approximation. The fragmentation is performed 
using the Lund string fragmentation scheme as implemented in JETSET [6] , which is also used 
to simulate the hadronisation process. The strangeness suppression factor in JETSET is chosen 
as 0.3, while a value asirriz) = 0.12 is used for the strong coupling constrant. 

From a study of Monte Carlo events, the hard component is found to be larger than the 
soft component. Their ratio goes from around two at low values of pt to around three at high 
values of Pf. 

^^If not mentioned otherwise, the symbol A refers to both A and A, as both charge conjugate final state are 
analysed. 
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The following Monte Carlo generators are used to simulate the background processes: 
KK2f [12] for the annihilation process e"'"e~ —>■ qq(7); KORALZ [13] for e"'"e^ r+r~(7); 
KORALW [14] for e+e" W+W" and DIAG36 [15] for e+e" e+e-r+r-. The response 
of the L3 detector is simulated using the GEANT [16] and GHEISHA [17] programs. Time 
dependent detector inefficiencies, as monitored during each data taking period, are included in 
the simulations. All simulated events are passed through the same reconstruction program as 
the data. 

3 Event selection 

Two-photon interaction events are mainly collected by the track triggers [18], with a low trans- 
verse momentum threshold of about 150 MeV, and the calorimetric energy trigger [19]. The 
selection of e"'"e~ —>■ e^e~hadrons events [20] is based on information from the central tracking 
detectors and the electromagnetic and hadronic calorimeters. It consists of: 

• A multiplicity cut. To select hadronic final states, at least six objects must be detected, 
where an object can be a track satisfying minimal quality requirements or an isolated 
calorimetric cluster of energy greater than 100 MeV. 

• Energy cuts. The total energy deposited in the calorimeters must be less than 40% of ^/s 
to suppress events from the e"'"e~ qq(7) and e+e" r+r^(7) processes. In addition, 
the total energy in the electromagnetic calorimeter is required to be greater than 500 MeV 
in order to suppress beam-gas and beam-wall interactions and less than 50 GeV to remove 
events from the annihilation process e^e~ qq(7)- 

• An anti-tag condition. Events with a cluster in the luminosity monitor, which covers the 
angular region 31 mrad < 9 <Q2 mrad, with an electromagnetic shower shape and energy 
greater than 30 GeV are excluded from the analysis. In addition, events with an electron 
scattered above 62 mrad are rejected by the energy cuts. 

• A mass cut. The invariant mass of all visible particles, Wms-, must be greater than 5 GeV. 
In this calculation, the pion mass is attributed to tracks while isolated electromagnetic 
clusters are treated as massless. The distribution of W^s for data and Monte Carlo 
after all other cuts are applied is shown in Figure [TJ Values of W^s up to 200 GeV are 
accessible. 

About 3 million hadronic events are selected by these criteria with an overall efficiency of 
45%. The background level of this sample is less than 1% and is mainly due to the e"'"e~ — >■ qq(7) 
and e"^e~ — > e^e~T^T~ processes. The background from beam-gas and beam-wall interactions 
is found to be negligible. 

The A baryons are reconstructed using the decay A pvr. Secondary decay vertices are 
selected which are formed by two oppositely charged tracks. The secondary vertices must satisfy 
the following criteria: 

• The distance, in the plane transverse to the beam direction, between the secondary vertex 
and the primary e"'"e~ interaction point must be greater than 3 mm. 

• The angle between the total transverse momentum vector of the two outgoing tracks 
and the direction in the transverse plane between the primary interaction point and the 
secondary vertex must be less than 100 mrad. 
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The distributions of these variables are presented in Figure |2l The proton is identified as 
the track with the largest momentum. Monte Carlo studies show that this association is correct 
for more than 99% of the vertices. In addition, the dE/dx measurement of both proton and 
pion candidates must be consistent with this hypothesis with a confidence level greater than 
0.01. 

After these cuts, about 70000 events are selected. The distribution of the invariant mass 
of the pvr system, m(p7r), shows a clear A peak, as shown in Figures El and ID for the different 
Pt bins listed in Table ^ The resolution of m(p7r) is found to be around 3 MeV and is well 
reproduced by Monte Carlo simulations. 

4 Differential cross sections 

The differential cross sections for A production as a function of pt and |?7| are measured for 
W^y > 5 GeV, with a mean value of 30 GeV, and a photon virtuality < 8 GeV^ with 
(Q2) ~ 0.2 GeVl This phase space is defined by cuts at the generator level of the Monte 
Carlo. 

The number of A baryons in each pt and |?7| bin is evaluated by means of a fit to the 
m(p7r) spectrum in the interval 1.085 GeV < m(p7r) < 1.145 GeV. The signal is modelled with 
a Gaussian and the background by a fourth-degree Chebyshev polynomial. All parameters, 
including the mass and width of the peak, are left free. The results are given in Tables ^ |21 
andlH 

The overall efficiencies, also listed in Tables ^ El and El include reconstruction and trigger 
efficiencies and take into account the 64% branching fraction of the decay A ^ pvr. The recon- 
struction efficiency, which includes effects of the acceptance and the selection cuts, is calculated 
with the PHOJET and PYTHIA Monte Carlo generators. As both generators reproduce well 
the shapes of the experimental distributions of hadronic two-photon production [20], their av- 
erage is used to calculate the selection efficiency. The efficiency does not depend on the 
cut-off. The track trigger efficiency is calculated for each data taking period by comparing the 
number of events accepted by the track and the calorimetric energy triggers. The efficiency 
of the higher level triggers is measured using prescaled events. The total trigger efficiency 
increases from 82% for pt < 0.4 GeV to 85% in the high pt region. 

The differential cross section as a function of pt is first measured for the three different data 
samples collected in 1998, 1999 and 2000 at different values of a/s and for different trigger and 
machine background conditions. Good agreement is obtained between the three measurements. 
The different data samples are therefore combined into a single data sample. 

Sources of systematic uncertainties on the cross section measurements are: background 
subtraction, scale and resolution uncertainties, Monte Carlo modelling and the accuracy of the 
trigger efficiency measurement. Their contributions are listed in Table IH The dominant source 
of systematic uncertainty is due to background subtraction. It is estimated using different back- 
ground parameterizations and fit intervals in the fitting procedure. The scale and resolution 
uncertainties are assessed by varying the selection cuts. The main contributions arises from 
the secondary vertex selection (3.2%) and the proton and pion identification criteria (2.5%). 
The uncertainty due to the selection of e^e~ e^e~hadrons events is 1%. The Monte Carlo 
modelling uncertainty, taken as half the relative difference between PHOJET and PYTHIA, 
increases with pt from 1% to 5%. A systematic uncertainty of 2% is assigned to the determi- 
nation of the trigger efficiency, which takes into account the determination procedure and time 
stability. 
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5 Results 



The average multiplicity of A baryons in the range W^^ > 5 GeV, 0.4 GeV < pt < 2.5 GeV and 
\rj\ < 1.2 is (1.57 ± 0.11) x 10~^ per e"'"e~ — > e^e~hadrons event. The uncertainty includes 
both statistical and systematic uncertainties. This result is below the value (1.80 ±0.01) x 10~^ 
predicted by PHOJET and above the value (1.43 ± 0.01) x 10"^ predicted by PYTHIA. The 
ratio of A to A baryons, as determined from the charge of the most energetic track of the vertex, 
is found to be N(A)/N(A) = 0.99 ± 0.04. 

The sum of the differential cross sections for the e"'"e~ — > e"''e~AX and e"^e~ — > e"*"e~AX 
processes as a function of pt for \r]\ < 1.2 is presented in Table ^ and Figure El Mass ef- 
fects explain the lower value obtained in the first bin. The behaviour of the cross section for 
0.75 GeV < Pt < 2.5 GeV is well described by an exponential of the form Aexp{—pt/{pt)), 
as shown in Figure with a mean value (pt) = 368 ± 17 MeV. For comparison, the values 
(pf) ~ 230 MeV and (pt) — 290 MeV are obtained for inclusive 7r° and Kg production, respec- 
tively [1]. 

The data are compared to the PHOJET and PYTHIA Monte Carlo predictions in Figure 
EJd. The region pt < 0.6 GeV is well described by PYTHIA, whereas PHOJET overestimates 
the cross section. On the contrary, the region 0.6 GeV < < 1 GeV is better reproduced by 
PHOJET while PYTHIA lies below the data. For pt > 1 GeV both PYTHIA and PHOJET 
underestimate the data. This level of agreement is not unusual in two-photon physics. 

The differential cross sections as a function of |?7| for 0.4 GeV < < 1 GeV and 1 GeV < 
Pt < 2.5 GeV are given in Tables El and El and shown in Figure IHl Both Monte Carlo programs 
describe well the almost uniform rj shape, while the size of the discrepancy on the absolute 
normalization depends on the pt range. 
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Pt 

(GeV) 


iPt) 
(GeV) 


Number of A 


Mass of A 
(MeV) 


Efficiency 

(%) 


da / dpt 
(pb/GeV) 


0.4-0.6 


0.50 


3412 ± 71 


1116.0 ±0.1 


10.2 ±0.1 


273 ± 6 ± 36 


0.6-0.8 


0.69 


4408 ± 85 


1115.9 ±0.1 


13.7 ±0.1 


264 ± 5 ± 29 


0.8-1.0 


0.89 


3420 ± 81 


1116.1 ±0.1 


15.3 ±0.2 


183 ± 4 ± 15 


1.0-1.3 


1.12 


3201 ± 87 


1115.7±0.1 


16.8 ±0.2 


104 ± 3 ± 7 


1.3-1.6 


1.43 


1222 ± 55 


1115.5 ±0.1 


17.7 ±0.4 


38 ± 2 ± 3 


1.6-2.0 


1.77 


578 ± 41 


1115.0 ±0.2 


15.9 ±0.6 


15 ± 1 ± 1 


2.0-2.5 


2.21 


292 ± 25 


1115.9 ±0.3 


17.2 ± 1.3 


6 ± 1 ± 1 



Table 1: The number of A baryons estimated by the fit, together with the A mass, the overall 
efficiency and the corresponding differential cross section as a function of for \r)\ < 1.2. The 
first uncertainty on the cross section is statistical and the second systematic. 



\v\ 


Number of A 


Efficiency 


da/d\rj\ 


0.0-0.3 


2953 ± 72 


13.8 ± 0.2 


59 ± 1 ± 3 


0.3-0.6 


2742 ± 63 


13.4 ± 0.2 


56 ± 1 ± 3 


0.6-0.9 


2904 ± 70 


13.0 ± 0.2 


61 ± 1 ± 4 


0.9-1.2 


2774 ± 89 


11.1 ± 0.1 


68 ± 2 ± 8 



Tabic 2: The number of A baryons estimated by the fit, together with the overall efficiency and 
the corresponding differential cross section as a function of pseudorapidity for 0.4 GeV < pt < 
1 GeV. The first uncertainty on the cross section is statistical and the second systematic. 



1^1 


Number of A 


Efficiency 
(%) 


da/d\ri\ 
(pb) 


0.0-0.3 


1458 ± 60 


18.8 ± 0.5 


21 ± 1 ± 1 


0.3-0.6 


1411 ± 62 


18.2 ± 0.4 


21 ± 1 ± 1 


0.6-0.9 


1480 ± 62 


19.3 ± 0.5 


21 ± 1 ± 1 


0.9-1.2 


1007 ± 63 


14.3 ± 0.4 


19 ± 1 ± 3 



Table 3: The number of A baryons estimated by the fit, together with the overall efficiency and 
the corresponding differential cross section as a function of pseudorapidity for 1 GeV < pt < 
2.5 GeV. The first uncertainty on the cross section is statistical and the second systematic. 
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Pt 


Background 


Scales and 


Monte Carlo 


Trigger 


Total 


(GeV) 


subtraction (%) 


resolutions (%) 


modelling (%) 


efficiency (%) 


(%) 


0.4-0.6 


12.1 


4.2 


1.4 


2.0 


13.0 


0.6-0.8 


9.8 


4.2 


1.6 


2.0 


11.0 


0.8-1.0 


6.6 


4.2 


1.9 


2.0 


8.3 


1.0-1.3 


5.3 


4.2 


2.1 


2.0 


7.3 


1.3-1.6 


5.8 


4.2 


2.4 


2.0 


7.8 


1.6-2.0 


6.4 


4.2 


2.8 


2.0 


8.4 


2.0-2.5 


19.4 


4.2 


4.8 


2.0 


20.5 



Table 4: Systematic uncertainty on the cross section of the e+e" — > e+e~AX and e'^e" — > 
e'*"e~AX processes due to background subtraction, scales and resolution uncertainties, Monte 
Carlo modelling and trigger efficiency. The total systematic uncertainty is the quadratic sum 
of the different contributions. 
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Data 

Pythia + bkg 








50 1 00 1 50 
W.i, (GeV) 



200 



Figure 1: The distribution of the visible mass W^is- The Monte Carlo distributions are nor- 
malized to the data luminosity. Various contributions to the background (bkg) are shown as 
cumulative histograms. 
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Figure 2: Distribution of the variables used for the selection of secondary vertices: a) the 
distance, in the plane transverse to the beam direction, between the secondary vertex and the 
primary e"'"e~ interaction point, d_L, and b) the angle between the total transverse momentum 
vector of the two outgoing tracks and the direction in the transverse plane between the primary 
interaction point and the secondary vertex, a. In each plot, all other selection criteria are 
applied. The predictions of the PHOJET Monte Carlo are shown as the full line and the 
contribution due to A baryons as the dashed line. The Monte Carlo distributions are normalized 
to the data luminosity. 
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Figure 3: The invariant mass of the pn system for a) 0.4 GeV < pt < 0.6 GeV, b) 0.6 GeV < 
Pt < 0.8 GeV, c) 0.8 GeV < pt < 1 GeV and d) 1 GeV < < 1.3 GeV . The signal is modelled 
with a Gaussian and the background by a fourth-degree Chebyshev polynomial. 
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Figure 4: The invariant mass of the pn system for a) 1.3 GeV < pt < 1-6 GeV, b) 1.6 GeV < 
Pt < 2 GeV and c) 2 GeV < pt < 2.5 GeV. The signal is modelled with a Gaussian and the 
background by a fourth order Chebyshev polynomial. 
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Figure 5: The differential cross section as a function of pt for the e+e~ — > e+e~AX and e+e~ — > 
e+e~AX processes for |?7| < 1.2: a) with the exponential fit described in the text and b) 
compared to the predictions of the PHOJET and PYTHIA Monte Carlo programs. 




Figure 6: The differential cross section as a function of |?7| for the e+e e~^e AX and 
e+e~ e+e^AX processes for: a) 0.4 GeV < pt < I GeV and b) 1 GeV < pt < 2.5 GeV. The 
data are compared to the predictions of the PHOJET and PYTHIA Monte Carlo programs. 
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